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^ • 1. Introduction 

pH I The aim of this paper is to review some of recent developments in mod- 

elling the Charge Current (CC) neutrino interactions with both free nucleons 
and nuclei targets in the few GeV neutrino energy region pp. This energy 
^ ' range is characteristic for atmospheric neutrinos and for several running or 

approved long-baseline experiments. The knowledge of the cross sections is 
necessary for future more precise measurements of neutrino oscillation pa- 
rameters. In our discussion we adopt a practical approach and will always 
have in mind Monte Carlo (MC) implementation of presented models. 

The few GeV energy region is rather complicated because three different 
dynamical formalisms are relevant: of quasi-elastic reactions, resonance ex- 
citation and more inelastic channels treated together in the DIS formalism. 
The significance of three dynamics is seen in Fig. 1 where the total CC cross 
sections for the muon neutrino scattering off free isoscalar nucleon target 
(i.e. average from proton and neutron targets)is presented. The contribu- 
tions from quasi-elastic, single pion production (SPP) and more inelastic 
channels (denoted as DIS) are also shown separately. It is seen that in the 
few GeV energy region all three contributions are important. The cross 
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Fig. 1. Total cross section for i/^ CC scattering on isoscalar target as predicted by 
the WROCLAW MC generator. The contributions from quasi-elastic, single pion 
production (SPP) and more inelastic channels (denoted as DIS) are also shown 
separately. 

section from SPP channels in Fig. 1 is restricted by the condition on the 
invariant hadronic mass W < Wcut = 2 GeV. 

The plan of the paper is the following. We start in Section 2 from a 
description of quasi-elastic reaction. In Sections 3 and 4 we review models 
of SPP and the formalism of Deep Inelastic Scattering. Sections 5 and 6 
deal with nuclear effects. 

2. Quasi-elastic reactions 

There are two CC AS* = quasi-elastic channels: ui + n ^ l~ + p and 
i^i + p ^ + n. In the discussed energy region the condition « My^ 
holds and it is enough to consider processes in the effective Fermi theory 
approximation . 

The matrix element contains leptonic part which is exactly known and 
the hadronic one, which cannot be calculated from first principles. The 
hadronic current contains four form- factors, functions of Q^. Vector form- 
factors Fi^2 are determined (CVC) by their electromagnetic counterparts. 
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Fig. 2. Quasi-elastic cross section for Ug,, Vfj, and Vr- Experimental points refer to 
v^j^ scattering. 

The form- factor Fp can be expressed in terms of Fa (PC AC) In recent 
years an improvement to MC codes was introduced by a replacement of old- 
fashioned dipole form-factors with one of the available fits to experimental 
data [H]. Fa is considered to be in the dipole form with two parameters: 
QA determined by the /? decay and the axial mass Ma which is not exactly 
known. The value of Ma determines the shape of da jdQ'^ and the overall 
quasi-elastic cross section 0]. 



3. Single pion production 

There are three CC SPP channels for neutrino reactions: f ; + n — > 
^~ + p + vr*^, i^i + n ^ ?~ + n + vr"'", z^; ^ Z~ + tt"^, and another three 
for anti-neutrino reactions. The characteristic feature of SPP reactions is 
the appearance of the A resonance in the differential cross section da jdW . 

There are several theoretical models of SPP. Almost all Monte Carlo 
generators use the Rein-Sehgal (RS) model It includes contributions 
from 18 resonances of mass M^es < 2 GeV treated in the coherent way. The 
non-resonant contribution is then added in the incoherent way in order to 
get an agreement with available data. An alternative model has been devel- 
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Fig. 3. Modification of the shape of da jdQ^ due to non-dipole electromagnetic form 
factors. For both curves Ma — 1-03 GeV. 
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Fig. 4. Modification of the shape of da/dQ^ due to different choices of the axial 
mass Ma — 1, 1.03, 1.07 GeV. It is seen that the choice of smaller Ma reduces the 
overall cross section in the different way than the substitution of dipole form-factors 
by BBBA05 ones (see the Fig. 3). 



oped recently by the Dortmund group [S]. It is based on experimental data 
on electromagnetic helicity amplitudes. The model contains the following 
resonances: P33(1232), Pii(1440), Di3(1520) and 511(1535). It includes mf 
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Fig. 5. Cross section for v^p fi pn^ as predicted by the WROCLAW MC 
generator. 



{mi is the charged lepton mass) terms which reduce the da/dQ"^ by ~ 5—10% 
at low and which are absent in the original RS model. A similar model 
with few resonances was constructed many years ago by Fogh-Nardulli 0. 

A common difficulty of SPP models is related to the issue of description 
of the non-resonant background. The most systematic approach is the one 
adopted in the Sato-Lee model 0. It is based on the quark model with the 
pion cloud effects taken into account. It predicts the non-resonant contri- 
bution to z^-neutron SPP channels on the level of ~ 25% and to z^-proton 
SPP on the level of ^ 5%. Alternative effective descriptions of non-resonant 
background are introduced in some MC codes. The idea is to simulate the 
background by a fraction of the DIS cross section in the resonance kine- 
matical domain i.e. for W < Wcut [S] ^01- Another theoretical possibility 
to deal with the non-resonant background is suggested by the hypothetical 
two-component quark-hadron duality jTJ. If the hypothesis is true, the 
background is given by the sea quark contribution to DIS structure func- 
tions. One can also try to model the non-resonant background by assuming 
that its dependence on the invariant mass is similar as in electron scattering 
i.e. ~ y/W^^^Wth^Yj<^j{Q'^){W - ^thrV , where Wthr is the threshold for 
the pion production |12j . 

The lack of precise experimental data makes it impossible to select a 
preferred model of SPP. 
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Fig. 6. Cross section for v^n — > ^ rnr^ as predicted by the WROCLAW MC 
generator. 




Fig. 7. Cross section for j/^n — > /i p7r° as predicted by the WROCLAW MC 
generator. 



4. More inelastic channels 

The DIS formalism describes the inclusive j^N cross section in the scaling 
limit. The justification of the theory comes from the perturbative QCD. In 
the few GeV neutrino energy region an important contribution to the cross 
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Fig. 8. Distribution of events in W and for Vfj_N scattering at energies E^, — 
1 GeV (top) and E^ = i GeV (bottom) as predicted by the WROCLAW MC 
generator. The quasi-elastic contribution is seen the peak at W = M (nucleon's 
mass). The A excitation region is also clearly seen. 



section comes from the small region, where the theory behind the DIS 
formalism is not valid. Thus the first problem is to get a correct form of 
the structure functions. The standard procedure is to express F^^^ in terms 
of Fi^2 !TB! and then to express Fi in terms of F2 and R = (cTL/aT)- In 
the scaling limit remaining i^2,3 are given as combinations of PDF's (parton 
distribution functions). In the region we focus on, the target mass and 
twist corrections must be taken into account 14 . The choice which is 
adopted in most MC codes is to apply structure functions with corrections 
modelled in analogy to the electron scattering case. Corrections which are 
available in the literature are applied to LO GRV98 PDF's ^Hj. Their form 
is closely related to the issue of quark-hadron duality: the DIS structure 
functions describe on average the electron scattering data in the resonance 
region 16 . It is an open problem whether quark-hadron duality should 
hold also in z/N scattering. There has been recently a lot of investigation 
in this field ^j. The theoretical analysis of the duality has been done in 
the framework of SU{6) quark model of resonances ^Sj. It is not clear 
whether arguments valid for the vector part of the hadronic current should 
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hold true also for its axial counterpart. Another problem is to understand 
what happens in the kinematical region < 0.5 GeV^, where one does not 
expect the quark-hadron duality to be present. In electron scattering it is 
known that -F|'^ — > and F^^ Q^, but for neutrino structure functions 
the presence of the axial current, which is not conserved makes the situation 
more complicated P^ . 

Once the inclusive cross section is calculated one has to evaluate contri- 
butions from exclusive channels. One possibility is to use the KNO theory 
which provides average multiplicities of particles in the final state ^Oj . The 
only remaining problem is then to redistribute to the particles energy and 
momentum transfer. Another strategy is to use the LUND fragmentation 
and hadronization routines j2()| I21j . 

In MC implementation of either scheme it is necessary to decide in sev- 
eral important points. Where should be a boundary between DIS and res- 
onance (SPP) contributions? The very definition of the RS model suggests 
that one should define Wcut = 2 GeV. Some authors argue that the RS 
model underestimates the cross section at higher W and the better choice 
is Wcut ~ 1-7 GeV ^2]- This is the choice adopted by the authors of 
the NEUGEN/GENIE MC code The comprehensive analysis of all 

the available date led other authors to the conclusion that one should take 
Wcut ~ 1-5 GeV This is approximately the choice implemented in 

the WROCLAW MC generator f^T where only A resonance contribution is 
included. 

5. Nuclear effects - generalities 

The treatment of nuclear effects depends on the neutrino energy. In 
the few GeV energy region one can rely on the picture in which neutrino 
interacts with individual nucleons inside nucleus (impulse approximation 
(lA)). It is not completely clear starting from which energies I A is the correct 
approach. Some authors argue that the Fermi gas model (the simplest 
mean field theory realization of lA) works well for electron neutrino energies 

> 200 MeV |23. Other authors are more conservative and argue that 
the lA picture makes sense for momentum transfers q > 400 MeV, which 
translates into higher neutrino energies [2^. For example, for neutrino 
energy Ei, = 0.8 GeV the contribution to the quasi-elastic cross section 
from the momentum transfers q < 400 MeV is ~ 20%. 

The simplest realization of the lA is known as PWIA (plane wave im- 
pulse approximation): one assumes that the nucleon produced in the pri- 
mary vertex leaves nucleus without further re-interactions. This is an obvi- 
ous oversimplification and it is better to include ESI (final state interactions) 
effects. A possible systematic approach to deal with FSI is known as DWIA 
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(distorted wave impulse approximation) [2SI- In MC codes FSI effects are 
usually treated by means of inter-nuclear cascade modules The propa- 
gation of nucleons, pions and other particles inside nucleus is semiclassical. 
It is important to implement the concept of the formation zone. Many other 
theoretical schemas to deal with FSI has been developed as well j27j. 

6. Nuclear effects - some models 

The advantage of the Fermi gas (FG) model is that it is easily applica- 
ble in MC routines. The basic FG model is defined by just two parameters: 
Fermi momentum and binding energy [^. It is necessary to deal with 
the problem of how to calculate off-shell nucleon matrix elements and the 
de Forest prescription is the common way to handle it [21] • One has also to 
decide about the kinematics. Smith-Moniz approach is the simplest choice 
and the other is to take into account the recoil nucleus momentum [HUj. An 
improvement to the above versions of the FG model is obtained in the frame- 
work of LDA (local density approximation): the Fermi momentum becomes 
a local quantity according to the density profile of the nucleus jHJ . Further 
improvement is introduced by modifying the momentum distribution of the 
nucleons by adding the high momentum tail JHH]. In the framework of the 
FG model the only FSI effect is Pauli blocking. 

The spectral function (SF) approach represents an improvement with 
respect to FG model by providing a realistic probability distribution of 
momenta and binding energies of nucleons inside nuclei j32j . Theoretical 
models of SF are obtained by combining the mean field (shell model) and the 
correlated part. The second one is relevant at higher values of momenta and 
binding energies and comes from short range interactions due to correlated 
high momentum pairs of nucleons j33j . Recently the correlated part of the 
SF has been directly measured in electron experiments |34j . 

Reliable models of SF exist for lighter nuclei up to oxygen The 
mean field part of SF is clearly seen as probability distribution peaks at 

E 42 MeV (Is), E 19 MeV (1^3/2) and E 13 MeV (1^1/2). The 

characteristic feature of the Is level is that it is smeared out in E. 

MC implementation of SF approach is straightforward. One can also 
use the effective approach in which the relevant information about the spec- 
tral function is contained in two functions: the probability distribution of 
nucleons momenta and the average momentum dependent binding energy 

IM|. 

In order to construct SF for heavier nuclei it is necessary to know the en- 
ergy levels and spectroscopic factors. The correlated part of SF is universal 
and depends only on the nucleus size [H^ . 
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Fig. 9. The effective description ^S] based on O. Bcnhar's oxygen spectral function 
|35| . The upper two plots show the momentum distribution (left-hand side) and 
the average momentum dependent binding energy (right-hand side). Below a com- 
parison is shown between three modelings of nuclear effects: by Fermi gas (FG), 
spectral function (SF) and the effective one (Eff). 

SF can be also used to model SPP in the resonance region. A new issue 
is the dependence of the A resonance width on the nuclear matter |37j . 

In the context of the DIS formalism there are specific methods to deal 
with nuclear effects to describe the shadowing, anti-shadowing etc. Recently 
a comprehensive model has been proposed to describe all the effects in the 
unique theoretical frame It includes also Fermi motion effects and in 
combining it with nuclear effects for quasi-elastic and SPP one has to be 
careful to avoid double counting. 

A step beyond lA would be to include contributions from 2-body current. 
Is is believed they are necessary in order to explain the excess of the cross 
section in the DIP region between quasi-elastic and A excitation peaks j39j . 
The problem with 2-body currents is that the computations which must 
be performed are algebraically very involved jlUI- Some authors tried to 
approximate the 2-body contribution to the neutrino cross section with the 
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conclusion that it is very important: for = 700 MeV it is responsible for 
up to ~ 25% of the cross section in the kinematical region of energy transfer 
uj G (80,250) MeV gU. 

7. Final remarks 

There has been efforts in the past to create a universal MC code to 
describe neutrino interactions So far each experiment uses its own MC 
focused on particular neutrino energy spectrum, target, detection techniques 
etc. The most promising ongoing project to construct a universal MC is that 
of GENIE 

All the theoretical considerations is this review were subject to big ex- 
perimental uncertainty in z^N cross sections. The good news are that the 
MINERvA experiment is under way [IH^. It will enable us to settle a lot of 
unknowns in free nucleon and nuclei targets cross sections. 
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